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Abstract 



The aim of this paper is to control the rate of convergence for central limit theorems of sojourn 
times of Gaussian fields in both cases; the fixed and the moving level. Our main tools are the 
Malliavin calculus and the Stein's method, developed by Nualart, Peccati and Nourdin. We also 
extend some results of Berman to the multidimensional case. 
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0< 

■ 1 Introduction 

Let X = {X{t),t £ R } be a stationary centered Gaussian field and T be a measurable subset of 
R"^. The sojourn time (or the volume of the excursion set) of X above the level ut in T is defined as 

/ KX{t) > UT)dt. 

J> . JT 

^ ■ The origin of this subject is the intersection between the study of the geometric properties of 

' random surfaces and the one of the non-linear functionals of Gaussian fields. Moreover, it has many 

04 ' applications in statistics of random processes (see, for example, Spodarev and Timmermann[B]). 

, The case of a fixed level: ut = u = const has been addressed in dimension 1 by of Sun |19) . 

I Chambers and Slud [3, Major and Giraitis and Surgailis [8]. Later on, some multidimensional 

versions were proved by Breuer and Major [J, Arcones pQ, Ivanov and Leonenko [5| and Bulinski, 
Spodarev and Timmermann 0. Their works are based on the following assumption 

' (A). {X{t) : t G R } is a stationary centered Gaussian field with unit variance and covariance 

5— ( ' function p(t) such that 

■ r 

■ ■ ■ / \p{t)\dt < oo, 

and can be presented in the following statement. 

Theorem 1. Let {X{t) : t £ R'*} be a random field satisfying the condition (A). For a fixed real- 
valued It, define the sojourn time as 



Then, as T tends to infinity, 



St = / MX{t) > u)dt. (1) 



■ AA(0,a 



where 



< = V ^ (")^"-^(") / p-it)dt < oo, (2) 



(y3 is the density function of the standard Gaussian law and $ is the tail of its distribution. 
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Berman [3] considered the problem for a Gaussian process in the case when the level depends 
on T. When the covariance function is not integrable, he assumed that the main component of the 
sojourn time is the first chaos in the Wiener chaos expansion. Else, his arguments were based on the 
spectral representation 

pW = / K't + s)b{s)ds, 

J R 

with the mixing condition b € L^nL^ and the m-dependent method. More precisely, he approximated 
the function & by a sequence of functions with compact support obtaining a family of m-dependent 
processes converges to the original one, and then he could use the central limit theorems that had 
been proved for this kind of process. His method can be applied in the multivariate case. 

However, the above works do not give us much information about the rate of convergence for the 
central limit theorems. Then, in this paper, we aim to control the speed in both cases: the fixed and 
the moving level. Our approaches come from the recent techniques, developed by Nualart, Peccati 
and Nourdin ([12|.[13|.|17|. etc.), that are the combination between the Malliavin calculus and the 
Stein's method. Here, we consider the Wasserstein distance for two integrable variables 

d{X,Y)= sup \E{hiX))-E(h(Y))\, 

heLip(l) 

where Lip(l) is the collection of all Lipschitz functions with Lipschitz constant < 1. Our main results 
are the following: 

Theorem 2 (Fixed level). Let {X{t) : t G R''} be a random field satisfying the condition (A). 
Assume that the covariance function p satisfies 

\p(t)\dt < (const) (log a)~ ^ , for a — >■ oo. (3) 



/M''\[-a,a]" 

Let St be defined by {7J). Then, 

d( ^"^^") ,AA(0,a^))<C(logr)-^^ 

where C is a constant depending on the field and the level, and satisfies 0). 
Note that the condition ((3} is weak, for example if 

p{t) ^ (const)\\t\\^°' , t — !> +00 

for some positive a > d , then it is met. Here and in the following, the notation f(x) = g{x), x — s> a 
fix) . 

means that lim ) : = 1. 

x-*a g(x) 

Theorem 3 (Moving level). Let {X(t) : t G M.''} be a random field satisfying the condition (A). 
Suppose that there exists a positive constant a £]0; 2] such that in a neighborhood o/O, the covariance 
function p satisfies 

1 - p(t) ^ (const)\\t\\°' for t -> 0. 
Let ut be a function that tends to infinity. One defines the sojourn time as 



St= KX{t) > UT)dt. 

Then, for every [3 G {Q;d/2), there exists a constant Cfj depending on the field such that 



\ (logr)i/6 ^ TP>^}(ut)ut 



In Nourdin et al |14) . the authors consider a very general case of Theorem [2] in the discrete time 
and obtain the bound under the form of an optimization problem. Here, in our particular case, we 
deal with a continuous time field and give an explicit bound. 
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2 Preliminaries 



In this paper, we use some notations that come from the Malliavin calculus introduced as follows. 

• Isonormal Gaussian process 

Let i5 be a real separable Hilbert space. Denote hy X — {X{h) : h G f)} an isonormal Gaussian 
process over Sj, that is a centered Gaussian family, defined on some probability space {Q, T , P), 
and E(X(/i)X((;)) = (ft, g)), for every h, g £ We assume that F is generated by X. 

• Wiener chaos expansion 

The n-th Hermite polynomial is 

Hr.{.)^{-^ire^j-,e^. 

For every n > 1, the n-th Wiener chaos Hn is defined as the closed linear subspace of -L^(f2, J^, P) 
generated by the random variables of the type H„{X{h)), where h £ is such that \\h\\s} — 1. 
Then, every square-integrable random variable Z G {Q, T , P) has the Wiener chaos expansion 

oo 

Z = ^J„(Z), (4) 

n=0 

where Jo{Z) = E(Z) and Jn{Z) is the projection of Z on Hn. Besides, for any n > 1 and 
h £ f), \\h\\sj = 1, the application 

J„(/i®") = H^{X{h)), 

can be extended to a linear isometry between the symmetric tensor product i^®" equipped with 
the norm \/nT|| .|| f,»„ and the n-th Wiener chaos %„. So, Z can be also decomposed in the form 



n = 



where /o(c) = c for all real c, /o = E(Z) and /„ £ , n > 1, are uniquely determined. 
Contraction and multiplication 

Let {efe, fc > 1} be a complete orthonormal system in S). Given / G i^®'' and g G i^®"*, then 
for every r = 0, 1, . . . ,p A g, the contraction of / and g of order r is the element of ^®(p+'?~2r) 
defined by 

oo 



Then, 



where f®rg G ^""^ is the symmetrization of / g),. g. 

Malliavin derivatives 

Let Z he a random variable of the smooth form 

Z = g{X{hr),...,XiK)), 

where n > 1, g : R" — >■ R is an infinitely differentiable function with compact support and 
hi G Then, the Malliavin derivative of Z is the element of L^(fl,^)) defined as 

i—l * 
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• Ornstein-Uhlenbeck operators 

The operator L is defined as i = X^J^o ~nJn- The domain of L is 

oo 

DomL = {Z e L'^{Q.) ■ ^n^\\J„{Z)\\l <oo}, 

n = l 

where ||Jn(-Z)||2 = || Jn(^)||L2(f2). Define the operator , caiied the pseudo-inverse of L, as 
L'^(Z) = V J„(Z) for all Z e L^iO). 

3 The fixed level case 

Lemma 1. For every n>2, let F„ be 

Fr. = -^ [ H4X{t))dt. 

Then, 

Var(||i3F„|||) < ^ J^^''-^' ^ ^) {2n - 2 - 2r)\ IpWIdf)' . 

Proof. The Malliavin derivative of F„ is 

DF,, = -^l nH„-j{X{t))i;tdt, 
VT'' J[o,T]d 

where is the element in Sj corresponding to X{t), i.e. X{t) — X{tl'(t)). And, 

\\DF„\\l = / s)i?„_i(X(t))J/„_i(X(s))dfds. 



[0,T]<'x[0,T]'' 



From the Mehler's formula, it is clear that 

l2i 1 2 



E[\\DF„\\i,]^—n' I {n-l)\p"{t-s)dtds = nVw[Fr,). 



Using the fact that 



[O,T]<ix[0,T]'' 



j/„_i(xw) = j„_i(V'f""^), 



and 



n-l 



7„_i(^f"-i)7„,i(^r-i) =Er!p /) /2„-2-2.(Vf"-'®.V'f"-') 



r = 
n-l 



2 



r = 



r 

2 



= E^' / p'-(t-^)/2n-2-2.Wf 



l-Dfrlli^ can be expressed as 

n-l 



,.=0 J[0,T]<ix[0,T] 

And, from the orthogonality of the chaos, the variance of ||_DFT||y-, is equal to 

n — 2 
T'^'^ " i[0,T] 

X (2n - 2 - 2r)!{V'f "''^'■(8)^'?""'^", V'®""'^''®'/'?""'"'') dtdsdt'ds' 
Each element of the scalar product has the form 

p"-'-^-'{t - t')p"-'-^-'is - s')p\t - s')p\s - t'), 



4 n — 2 „ 

/ {r\f{-~^Yp^^\t-s)p^+\t'^s') 



for some i £ [0; n — 1 — r]. And 

/ p'+\t - s)p''+\t' - s)p"~^~'"\t - t')p"-^-''-'{s - s')p'{t - (s - t') 

is at most equal to 

\p{t — s)p{t' — s')p{t — t')p{s — s')\ dtdsdt' ds' , 



[0,T] 



dx4 



/ \p{t — s)p{t' — s')p{t — s')p{s — t')\ dtdsdt' ds' 

Jr().Tl<ix4 



' [(),T]<* 

With the change of variable y = {t — s,t' — s' ,t — t' , s') 



\p{t — s)p{t' — s')p{t — t')p{s — s')\ dtdsdt'ds' 



dyi / \piyi)piy2)piy3)piy2 + ya - yi)\ d?/idj/2dj/3, 



can be written as 



where Ay^ is some domain in that depends on j/4. It is at most equal to 

/ dyj' \p{yi)p(y2)p{y3)\dy^dy2dy3 = T^ ( f \p{t)\dt). 

The same bound is obtained for the others. So, the variance of HDFrlly-, is at most equal to 
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□ 



In this paper, we use some facts about Hermite polynomials (see Szego \20\). 
Lemma 2. • For a fixed point u, there exists a constant Cu such that 

e-^'/^IH^Cu)! < C„(n/e)"/' V n e N. (5) 



There exists a constant K such that, for all u, n, 

ip{u)\Hn{u)\ 



As n tends to infinity, 



< K. (6) 



max e 



'^\H,,{x)\ ^ {const)V^.n-^^^^. (7) 



Proof of Theorem It is clear that 

+d [ '^■"--^^^"-\ U{0, a%^)) + d (AA(0, al^),m, a^)) ^ d, + d^ + d,, 

where St.n^ is the truncation of St at position A'^r in the Wiener chaos expansion. Nt will be chosen 
later on. 
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i) (Bound for d±) It is easy to show that 

f St — E(S't) St,Nt ^ E(St,jvt; 



d 



St — St,Nt 



/yd 



E 



Here, from (O and the Stirhng formula 

n! ~ \/27m(n/e)" 

we obtain the bound for di 



p"(t) Il (T - \t,\)dt. 



di < Cu^M^J [ \p{t)\dt 



J2 'i'^'^*' < {const)N:^'^\ 



(8) 



ii) (Bound for 0(2) From Theorem 3.1 of [12| . it is clear that 



/ St,Nt ~ E(St,jv7,) 



< 



CaTt, — ( ;= . — l^-L 



Nt 

p,5=l 



St - £(5-7) 



St - E{St] 



where Jp is the component in thep-th chaos defined in Q and (Jt.jj the variance of Jp[^^^—^£^] 
- Ifp = g= 1, 



o"T,i — ( DJj 



St — E(St) \ _ J / 5t — E(St 



Ifp=g> 1, 



2 -1 / m f — E(S't) 



St ~ E(5't 



,DJ„ 



^ p 



1^2 



Var 



T'* Jo 



Hp{X{t))dt 



Then, from Lemma [T] it is at most equal to 



If p > 1 and q — 1, then 



a /p-2 



p{t)\dt) (^(r!)2(P M (2p-2-2r)! 



T"^ i[0,T]<i 



Hi{X{t))dt 



Hp(X{t))dt ,D 



T'^ JlO,T]d 



P 

J^d 



[0,T]''x[O,T]<i 



p(t - s)Hp-i{X{t))dtds. 
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So, its variance is 



1 2 
P 



J'2d 



{p- iy.p(t- S)p{t' - s')p''-\t-t')dtdsdt'ds' < ^^^^ ^ 



'[0,T]''x 

If p, g > 1 and p q, then 
D ' 



W)\dt 



yd 



J'd / J 



p{t - s)Hp-i(X{t))Hg-i{X{s))dtds 

-l\fq-l 



[0,T]<'x[0,T]<' 
pAq — 1 



-'[o.r]''x[o,T]'' 
So, its variance is at most equal to 



p{t - s)Ip+,-2-2r{{^f^-^®r^f''-^)s)dtds. 



< 



{piY 



\p{t)\dt) Yl 



r = 



2 / \ 2 

n Iq-l 



(p + g-2-2r)! 



We obtain the bound for d2 



\p{t)\dt 



p=2 



(p!)2 yr^\J^ 



p-2 



P- 1 



(2p-2-2r)! 



+ E 

So, 



1^1^ ip\u)\Hp-^{u)H^-,{u}\ pq 

p q 



' p - 1 \ / g - 1 



d2 < [const) 



yd 



(p + g-2-2r)! 



(9) 



Indeed, from 



p-2 



((p-l)!)2 

p-2 

E 

r = 
p-2 

E 

r = 



r=0 
2 



WE('-Or;M (2P-2-20! 



p-l\ / 2p - 2 - 2r 
r / \ p — 1 — r 

2 

P ~ 1 ^ 22p-2-2r 

p-2 



< 22p-M^r ^12- 



p-i 



and the first term is at most equal to 



1 !X 3"-! 
(const)— }J : 



/yd^ P 



and the same for the second term, 
iii) (Bound for da) It is easy to show that 

d2(AA(0,air^),AA(0,cr2)) < {const){a'' - aj,^) 



Nt 

E 



3 = 1 



yd 



-dt 



^\u)Hl_^{u) 



+ V 

./H''\[-T,T]'' 



+ E 



„2/„,\ tt2 



I 



p''{t)dt 



n. 



p"{t)dt 
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From part i), the third term is at most equal to (const) Nj,^ ^'^ . For the first term, it is equal to 



d 



+ / , ^ w 

" i[-T,T]<i\[-VT,v^]'' 



71=1 



which is at most equal to 



^/^ 



( COTlst I 

The first part is at most equal to — -=h-. The sum of the second part and the second term is 



and at most equal to (const) (log T)^^ (from Q). So, 

da < {const){N-^^^ + T'^^'' + (logT)-^/^). (10) 
Summing up three bounds and H10|) . by choosing Nt = (logr)/4, we have the result. □ 

4 The moving level case 

In this section, we assume that the level depends on T and we denote by ut- Then the sojourn 
time 

St= f HX{t) > UT)dt 

i[0,T]<i 

has 

E(5't) = T'^'^{ut) 

^'^'^ (t) 

Vm{ST)= I n{T-\tj\)dtf ip{uT,UT,y)dy, 



ip(UT,UT, y) = 7== exp 



where 

27r\/l-y2 VI + y 

is the density of the bivariate normal vector 

' ' ^ y 
y 1 



When Ut tends to infinity. 



Var(gT) „ 
yd ^ 

then the Theorem [T] and [5] no longer hold. So, at first, we generalize the results of Berman [3] (chapter 
8) to estimate the variance of St (the detailed proofs are given in the Appendix). 
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Lemma 3. If the covariance function p satisfies the conditions in Theorem\^ then, for every e > 0, 
/ Ii{T-\tj\)dt •4){uT,UT,y)dy^T / / if{uT,UT,y)dydt, 

J[-T,T]d^-^ Jo i[-e,e]'' io 

for T, ut ^ oo. 

So, let B{u) be some function that satisfies 

fp(t) 



Then, 



[-e,E]<i Jo 



B{u) = I I (p{u,u,y)dydt, for it — >■ cxj. 

Var(S'T) ~ T'^B{ut), 



for T, ut — oo. 

Lemma 4. If the covariance function p satisfies the conditions in Theorem\^ then, 

B{u) = (const) ^^2+1 ' ^'^^ u ^ oo. 

Proof of Theorem [5l The distance between — - and the standard Gaussian variable is at 



most equal to 



VVar(ST) 



St — E(S't) St,Nt ~ E(5't,jvt) \ , j I St.Nt ^ E{St,Nt) 



\/Var(5T) ' \/Var(ST 



Y/Var(5T,]VT) 



,AA(o,i; 



where St.Nt is the truncate variable of St at position A^t in the Wiener chaos expansion. Nt will 
be chosen later on. 

• The first term is at most equal to (up to some multiplicative constants) 



Var(S'T - Snt) 



Var(S'T 
= (const) 



E„=iv^+i -\ J[-T,T]<i p - 

Var(5'T) 



= (const) 

< (const) 

< (const) 



E 

n = N-r + l 



v{ut)HI_^(ut) 



[-T,T]'i 



.n (r-li.l) 



\p(t)\dt 



ip(uT)Hl_^(uT) 



\ 



J2 n-^'^i 



n = Nrr + l 

where in the third line, we use the approximation 



Var(5'T) = T''B(ut) ^ (const)T''^^^^ , 



and in the last one, the fact Q is used. Then, we have the bound 
(const) 



B(ut) 



< (const) t 



(11) 
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For the second term, as the same argument in part ii) in the proof of Theorem [21 we have the 
bound 

(const) - 



which is at most equal to 

T'i/^lfi{uT)uT' ^^^^ 

Summing up ^ and by choosing Nt such that 3^"^ = y-/3+d/2^ ^j^g ^.gg^j^- follows. □ 

We have the following corollary 

Corollary 1. Let {X{t) : t G R''} be a random field satisfying the condition (A). Suppose that there 
exists a positive constant a g]0; 2] such that in a neighborhood o/O, the covariance function p satisfies 

1 - p{t) = (const) ||f II" for t 0. 

One defines the sojourn time 



St = / MX{t) > UT)dt. 

J[0,T]rf 

Let ut be a function that tends to infinity. Then, if 



(logT)-^/^^°° ^0, 



one has 

St - E(gT) d 
^/VMSt) ^ 



N{Q,1). 



2+Q 

Proof. Since (log T)"^/^iij,° — >■ 0, it is easy to see that 



Tl^ip{uT)uT ' 

for all P G (0; d/2) .From Theorem O the result follows. □ 

This extends, under the stronger hypothesis on ut, the results of Berman to Gaussian fields in 
E"* with d > 1. 



Appendix: Proofs of the Lemmas 3-4 

In this Appendix, we prove the Lemmas 3-4 analogously to the similar ones in with some minor 
changes. 

Proof of Lemma 3. It suffices to show that 

/[_TT1^\[-..1<< .n (T- lt,l)/;«0(^T,UT,J/)dydt 

'■ (13) 

n^(r - U)!^''' H^T,UT, y)dydt 

tends to for ut ,T ^ oo. In fact, denote 

77 = 1 - max(lp(s)| : s ^ {-e,tf). 
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If 77 = then there exists x such that \p{x)\ — 1, then the field is x- or 2x- periodic and the 
integral / |/3(t)|cit can not converge. Therefore, 77 is strictly positive. Since the function (/^(ut, iir, y) 
is increasing with respect to y, the numerator in (|13p is at most equal to 

tM^t,ut,1^v) f \p{t)\dt. (14) 

i[-T,T]'i\[-£,€]<i 



The denominator in (|13|1 can be decomposed as 

f n(T-|t,i)r" if{uT..UT,y)dydt- I Yl(T-\t,\)l .fi{uT,UT,y)dydt. (15) 

J[_E,e]dJ = l Jo J[_£,£]<13 = 1 J-p(t)- 



There exists a positive constant c < 1, such that p{t) < c, Vt, then the second term in H15|l is at 
most equal to 

(2e)'*T^-^=L=^^(nT). (16) 



Choose 5 < Tj and e' < e such that 



min(p(t) : t £ [-e ,eX) > 1 - 5, 



then the first term in (jlSp is lower-bounded by 

rp(t) 



f n (T-\tj\) f ip{uT,UT,y)dydt 

J[-C',c']d J = ^ Jo 
Jr_E' e'ld Jl- 



(■P{t) 

ip{uT,UT,y)dydt. 



and it has the lower bound 



{T -e'fip{uT,UT,l-S) [ {p{t) - 1 + 5)dt. (17) 

It is clear that p4|l and (|16p are negligible with respect to (|17p when ut and T tend to infinity, it 
implies the result. □ 

To prove the lemma |31 we need the following two results: 
Lemma 5. For every 6 > 1, there exists a constant K{6) > 0, such that, asymptotically 

B{u) > K{e)e^p{-ue/2). 

Proof. It suffices to prove the lemma for 9 in a neighborhood of 1. In such case, using p7p . we can 
choose 5 such that 

exp(— ?i^^/2) = ip{u, It, 1 — 5), 
and we are done. □ 

Lemma 6. For every 5 G (0, 1), one has 



B(u) 

limsup ^- — ^ < 1, 



2 



dt 



and 



liminf , ^ > 1. 



U—¥OC 



[2(2-5)]V2^r $(Jiz£mr/^)di 
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Proof. For any 5 e (0, 1), there exists e > such that 1 — p{s) < 5, Vs £ [— e, e]''. Then, 

/ V>{u,u,y)dydt ^ {2eY <p{u,u,y)dy < {2e)'^{l - 5)<p{u,u,l - 5). 

Jo Jo 



Since 



(^(u, u,l — 5) 



: exp 



2-S 



2^^1- (1-5)2 

and from Lemma[Sl B{u) is asymptotically greater than K{d) exp{—u'^9/2) for every 9 > 1, then by 
choosing 

i<e < 



2-5' 

/[_£ sjd /o^ * w, y)dydt is negligible with respect to -B(u) when u tends to infinity. Hence, B{u) is 
asymptotically equal to 

//■pW 
/ ip{u,u,y)dydt 
_E,£]£i Jl-S 



P(t) 



which is equal to, by the change of variable z — u'^{l — y), 



1-y 



l + y 



1/2 



dydt, 



V5 



2 - 2/1^2 



1/2 



dzdt. 



An upper bound of H18|l is 



^/2~~^ u 

which is equal to, by the change of variable x = \fzj2 

2^/2 ip{u) 



, dz 

V9(\A/2)^dt, 



2V2 ip{u) 



•(i-p(t)) 



(p{x)dxdt 



A lower bound of (fTg]) 



\/2 — 5 It 

V5(m) 



1 - p(f) 



11/2 



^ f 



2-5 



1/2 



dz 



2z 



dt, 



which is equal to, by the change of variable x = ^ z/{2 — 5), 



v/2(2 - 5) 



ii{«/{2-«))l/2 



ip{x)dxdt 



1 - 



-,1/2 



2-5 



<1> u 



1 1/2 



2-5 



dt. 



Since 



lim sup 

the lower bound is asymptotically equal to 
[2(2-5)]^/^^ 



■p(t) 



1/2 



1/2 



/ $ ( « 






2-(5 



0, 



-I 1/2 



dt. 



(18) 



□ 
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Proof of Lemma^ By change of variable t = z/u'^^°', the asymptotically upper bound in Lemma 
[6] is equal to 

It is clear that 

11^(1 — p{z/u^^'^)) — > C||2:||" for u — > oo, 
then by dominated convergence, this upper bound is asymptotically equal to 

By the same argument, the lower one in Lemma [6] is asymptotically equal to 

[2(2-^)1^/^41- / HC\\z\ndz. 

Let 5 tend to 0, we obtain the result. □ 
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